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Knowledge  of  atmospheric  turbidity  coefficients  is  very  important  in  meteorology,  climatology, 
atmospheric  pollution  monitoring,  and  in  the  prediction  of  solar  energy  availability  under  cloudless 
skies.  This  paper  provides  a  thorough  review  on  a  number  of  atmospheric  turbidity  indices  and  on  the 
several  methods  that  have  been  developed  in  the  past  few  decades.  The  Lambert-Bouguer-Beer  law  is 
the  basic  relationship  underlying  the  derivation  of  various  turbidity  indices.  Turbidity  can  be  referred  to 
as  monochromatic  (narrow  band)  wavelengths,  broadband  wavelengths,  and  the  total  spectrum. 
Narrow  band  turbidity  is  measured  using  sun  photometers,  while  the  other  two  turbidities  are 
measured  using  pyrheliometers  with  broadband  pass  filters.  The  Angstrom's  turbidity  coefficients 
and  the  Linke’s  turbidity  factor  are  among  the  most  frequently  used  atmospheric  turbidity  coefficients. 
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1.  Introduction 

Most  energy  produced  in  the  fusion  furnace  of  the  sun  is 
radially  transmitted  as  electromagnetic  radiation,  which  is  popu¬ 
larly  called  sunshine  or  solar  energy  [  1  ].  The  sun  has  a  tempera¬ 
ture  of  5760  K,  and  it  emits  the  total  power  of  approximately 
3.8  x  1030 *  W  or  6.25  x  1011  W  m  2.  Since  the  distance  between 
the  sun  and  the  earth  is  large,  the  amount  of  solar  radiation  that 
reaches  the  outside  of  the  earth’s  atmosphere  is  quite  low  at  only 
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1367  W  m~2,  which  is  called  the  total  solar  irradiance  (TSI)  or  the 
solar  constant  [2],  The  solar  constant  is  measured  on  the  surface 
perpendicular  to  the  sun’s  rays  at  the  average  sun-earth  distance 
above  the  atmosphere  [3].  The  intensity  of  solar  radiation  above 
the  atmosphere  approximately  varies  sinusoidally  over  the  year 
with  amplitude  close  to  3.3%  of  the  solar  constant  and  it  achieves 
its  maximum  around  the  first  day  of  January.  This  variation  arises 
from  the  variation  in  the  distance  of  the  earth  from  the  sun  [4]. 
Variation  in  solar  radiation  is  also  caused  by  the  difference  in  the 
emission  intensity  from  the  sun  itself  [3], 

Not  all  solar  radiation  received  at  the  periphery  of  the  atmo¬ 
sphere  reaches  the  Earth's  surface  [5].  Solar  radiation  is  partially 
depleted  and  attenuated  as  it  traverses  the  atmospheric  layers  [6]. 
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Nomenclature 

SD(z) 

the  relative  optical  thickness  associated  with  aerosol 

extinction  and  gaseous  absorption  other  than  ozone 

P 

the  Angstrom’s  turbidity  coefficient 

SRK(mr) 

the  integral  Rayleigh  optical  thickness  computed  by 

CL 

the  Angstrom  wavelength  exponent 

the  formula  of  Kasten  [63] 

spectral  aerosol  optical  depth  at  the  wavelength  X 

Tudmr ) 

Linke’s  turbidity  factor  determined  by  Kasten's  algo- 

Vis 

visibility  in  kilometers 

rithm  [63] 

Id 

diffuse  solar  radiation  (W  m~2) 

TLa(mr ) 

Linke’s  turbidity  factor  determined  by  Kasten's  algo- 

I„ 

direct  solar  radiation  (W  m~2) 

rithm  [68] 

h 

the  irradiance  of  the  extraterrestrial  solar  radiation 

SRa(mr) 

integral  Rayleigh  optical  thickness  computed  by  the 

he 

extraterrestrial  solar  irradiance  at  normal  incidence 

formula  of  [68] 

(1367  Wm-2) 

Elvis 

Linke’s  turbidity  factor  in  the  visible  region  of  the 

Ics 

threshold  irradiance 

solar  spectrum 

I  bn 

normal  direct  solar  radiation  at  the  earth’s  surface 

Ilm 

direct  hourly  normal  solar  radiation 

(Witt2) 

Tup 

Linke  turbidity  factor  obtained  using  Pinazo  et  al. 

Eo 

correction  term  for  the  earth-sun  distance 

method  (dimensionless) 

I*bn 

normal  incidence  solar  radiation  at  the  bottom  of  a 

Tr 

the  mean  Rayleigh  extinction 

dust-free  atmosphere 

TW 

the  selective  absorption  coefficient  for  atmospheric 

T5 

the  aerosol  optical  thickness  at  0.5  pm 

gases 

Pg 

Angstrom  turbidity  coefficient  obtained  using  Guey- 

Tu 

Unsworth-Monteith  turbidity  coefficient 

mard  method  (dimensionless). 

T[AM2 

Linke’s  turbidity  factor  for  an  air  mass  equal  to  2 

mr 

relative  optical  air  mass  at  a  given  pressure 

w 

precipitable  water  thickness  (cm). 

(dimensionless) 

I<db 

the  standardised  value  (to  zero  altitude  and  total 

Tpin 

Linke  turbidity  factor  obtained  using  Pinazo  et  al. 

column  ozone  equal  to  0.3434  atm  cm)  of  the  ratio 

method  (dimensionless) 

between  diffuse  and  direct  irradiance 

Kdb 

the  ratio  between  the  diffuse  and  direct  beam  normal 

B 

Schiiepp’s  turbidity  coefficient 

irradiance 

m* 

an  absolute  air  mass 

4aff 

the  effective  wavelength  corresponding  to  each  of  the 

CDA 

the  number  of  clean-dry  atmospheres 

filters 

7 

solar  elevation  angle 

the  aerosol  optical  thickness  at  a  specified 

Tu 

illuminance’s  turbidity  factor 

spectral  band 

Tv 

Volz’s  turbidity  factor 

Tl 

Linke’s  turbidity  factor  (dimensionless) 

*5 

Rayleigh  scattering  coefficient  for  air  moleculares  at 

Elm 

monthly  average  values  of  Linke’s  turbidity  factor 

the  wavelength  X 

SR(mr) 

the  integral  Rayleigh  optical  thickness 

Toz 

absorption  coefficient  for  ozone  at  X 

<5  r(z) 

the  relative  optical  thickness  relating  to  Rayleigh 

Tk 

Kastrov’s  turbidity  coefficient 

scattering  by  the  gaseous  molecules 

Tr 

rational’s  turbidity  factor 

The  presence  of  aerosol  particles  in  the  atmosphere  is  one  of  the 
most  important  factors  that  affect  the  amount  of  solar  radiation 
that  reaches  the  Earth’s  surface  under  cloudless  skies  [7].  Aerosols 
are  small  particles  that  are  either  in  the  solid  or  liquid  state  [8]. 
Common  aerosol  particles  have  radius  ranging  from  1  to  105  nm, 
while  very  small  particles  (called  Aitken  particles)  range  from  1  to 
10  nm  and  large  particles  range  from  102  to  103  nm.  Particles  in 
the  range  of  103  to  105  nm  are  called  giant  particles  [9,10].  Based 
on  the  formation  of  aerosol  dispersions,  aerosol  particles  have 
two  categories,  namely,  dispersed  and  condensed  aerosols  [8]. 
Aerosols  have  two  dominant  layers  in  the  atmosphere.  The  first 
dominant  layer  is  near  the  Earth’s  surface  at  0-3  km,  which  is 
affected  by  natural  dust  storms  and  man-made  inputs  to  the 
atmosphere,  whereas  the  second  dominant  layer  is  the  strato¬ 
spheric  dust  layer  at  15-25  km  above  sea  level,  which  is  caused 
by  volcanic  action  and  cosmic  sources  [11]. 

Aerosols  have  direct  radiative  forcing  because  they  scatter  and 
absorb  solar  and  infrared  radiations  in  the  atmosphere.  Aerosols 
also  alter  warm,  ice,  and  mixed-phase  cloud  formation  processes 
by  increasing  droplet  concentrations  and  ice  particle  concentra¬ 
tions,  consequently  causing  indirect  radiative  forcing  associated 
with  the  changes  in  the  microphysical  and  optical  properties  of 
clouds.  Furthermore,  aerosols  are  highly  absorbent  of  solar  radia¬ 
tion,  such  as  black  carbon,  and  they  may  reduce  cloud  cover  and 
liquid  water  contents  by  heating  the  cloud  and  the  environment 
where  the  cloud  forms.  This  is  known  as  the  semi-direct  effect 
because  it  is  the  result  of  the  direct  interaction  of  aerosols  with 


radiation  that  indirectly  influences  the  climate  by  altering  clouds. 
The  net  effect  of  aerosols  on  climate  depends  on  the  summation 
of  all  the  three  mechanisms  above  [12].  An  atmosphere  that 
contains  aerosols  is  described  as  turbid  or  hazy.  An  aerosol-laden 
atmosphere  has  the  property  called  atmospheric  turbidity,  which 
depletes  direct  solar  radiation  [10]. 

Studies  on  atmospheric  turbidity  are  important  in  meteorology, 
climatology,  and  in  atmospheric  pollution  monitoring.  Atmospheric 
turbidity  parameters  are  necessary  for  determining  the  amount  of 
spectral  global  irradiance  and  for  designing  photovoltaic  cells  and 
selective  absorbers  for  spectral  thermal  collectors  [13].  Aerosols 
make  the  atmosphere  turbid,  and  thus,  the  presence  of  aerosols  is 
commonly  quantified  by  an  index  of  turbidity.  In  the  last  century, 
some  indices  of  turbidity  were  proposed  and  several  techniques 
were  developed  to  determine  their  values  [14].  Among  these  indices 
are  the  Angstrom’s  and  the  Schiiepp’s  turbidity  coefficients,  which 
are  spectral  in  nature,  and  the  Unsworth-Monteith’s  coefficient 
and  Link’s  turbidity  factor,  which  are  broadband  in  nature  [15]. 
This  study  mainly  aims  to  investigate  atmospheric  turbidity 
and  then  describe  available  techniques  for  calculating  various 
turbidity  indices. 


2.  Atmospheric  aerosol  characteristics 

The  word  aerosol  was  introduced  more  than  80  years  ago  as  an 
analogy  to  the  term  hydrosol,  a  stable  liquid  suspension  of  solid 
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particles.  Aerosols  are  products  of  a  complicated  totality  of 
chemical  and  physical  processes  [16,17].  Aerosols  are  defined  in 
their  simplest  forms  as  small  solid  or  liquid  particles  that  remain 
suspended  in  the  air  and  follow  the  motion  of  the  air  within 
certain  broad  limits.  An  atmosphere  that  contains  aerosols  is 
described  as  turbid  or  hazy.  An  aerosol-laden  atmosphere  has  the 
property  called  atmospheric  turbidity,  which  depletes  direct  solar 
radiation  [16,10],  The  presence  of  aerosols  in  the  Earth’s  atmo¬ 
sphere  has  important  effects  on  the  transmission  of  solar  radia¬ 
tion  and  on  the  radiative  heat  transfers  in  the  atmosphere.  This  is 
due  to  the  fact  that  aerosols  absorb  and  scatter  solar  radiation  as 
it  passes  through  the  atmosphere  [14].  Scattering  refers  to 
processes  in  which  photons  change  direction  after  an  interaction, 
while  absorption  occurs  when  photons  are  removed  from  a  beam 
light  and  their  energy  is  converted  to  an  excitation  of  atoms  or 
molecules  [18].  Aerosols  are  physically  divided  into  the  following 
three  categories  according  to  particle  size: 

1.  Aitken  nucleus  with  particle  radius  smaller  than  101 2  nm 

2.  Large  nucleus  with  particle  radius  larger  than  102nm  and 
smaller  than  103 4  nm 

3.  Giant  nucleus  with  particle  radius  larger  than  103  nm 

Based  on  the  formation  of  the  aerosol  dispersions,  aerosols  fall 
into  the  following  two  categories  [8[: 

1.  Dispersed  aerosols:  the  aerosol  materials  that  are  either  in  the 
solid  or  liquid  state  are  pulverized  into  a  granular  state  via 
mechanical  pulverization  or  natural  weathering,  and  are  then 
suspended  in  the  air  because  of  the  increase  in  wind  force. 

2.  Condensed  aerosols:  the  formation  of  the  aerosol  particle  or 
droplet  is  caused  by  the  condensation  of  super-saturated  gas 
on  the  condensation  nucleus  or  by  the  mixture  of  different 
gases  through  photochemical  reaction. 

Studies  on  the  physical  properties  of  aerosol  particles  are 
important  for  the  assessment  of  the  effects  aerosol  particles  on 
the  climate  and  for  the  development  of  more  accurate  remote¬ 
sensing  procedures  from  satellite  sensors  [19].  Aerosol  properties 
vary  with  time  and  from  one  region  to  another  because  of  their 
short  lifespan  and  highly  temporal  and  spatial  variability  mea¬ 
surements  [20,21].  The  spatial-temporal  variability  of  aerosol 
characteristics  is  large  and  the  observational  data  are  so  frag¬ 
mentary  that  the  total  content  of  various  types  of  aerosols  is 
difficult  to  estimate,  while  the  available  power  estimates  of  the 
global  sources  of  natural  and  anthropogenic  aerosols  are  rather 
approximate  [17]. 

Kondratyev  et  al.  [17]  and  Giannakaki  et  al.  [22]  reported  that 
the  types  of  aerosols  depend  on  their  sources  and  are  roughly 
categorized  into  anthropogenic  or  natural.  Anthropogenic  sources 
are  determined  by  human  activities,  such  as  industrial  wastes 
from  chimneys,  toxic  exhausts  from  cars,  fires,  explosions,  soil 
erosion  in  agriculture,  and  open  mining.  The  estimates  of  anthro¬ 
pogenic  aerosols  are  more  reliable  compared  with  those  of 
natural  aerosols,  especially  in  regions  along  the  World  Ocean 
and  in  the  continents  that  are  difficult  to  access.  Atmospheric 
aerosols  are  classified  into  the  following  types  of  natural  aerosols 
based  on  their  composition  or  sources: 

1.  Products  of  sea  spray  evaporation. 

2.  Mineral  dust  wind-driven  to  the  atmosphere. 

3.  Volcanic  aerosols  (both  directly  emitted  to  the  atmosphere 
and  formed  by  gas-to-particle  conversion). 

4.  Particles  of  biogenic  origin  (directly  emitted  to  the  atmosphere 
and  formed  by  the  condensation  of  volatile  organic  compounds, 


such  as  terpenes,  and  from  the  chemical  reactions  between 
these  compounds). 

5.  Smokes  from  biota  burning  on  land. 

6.  Products  of  natural  gas-to-particle  conversion  (e.g.,  sulfates 
from  reduced  sulfur  on  the  ocean  surface  with  emissions  of 
dimethyl  sulfide). 

Important  types  of  anthropogenic  aerosols  also  include  the 
following: 

1.  Direct  industrial  emissions  of  particles  (soot,  smoke,  and  road 
dust,  among  others). 

2.  Products  of  gas-to-particle  conversion. 

The  amount  of  aerosols  in  the  atmosphere  is  sometimes 
measured  in  terms  of  the  number  of  particles  per  cubic  centi¬ 
meter.  The  number  of  dust  particles  is  usually  greater  on  land 
than  over  water  and  during  the  drier  seasons.  Moreover,  the 
number  of  dust  particles  is  lower  in  cold  polar  or  arctic  air  and  is 
higher  in  tropical  air.  The  number  density  exponentially  decreases 
with  altitudes  of  up  to  approximately  5  km  and  remains  constant 
at  altitudes  of  10-15  km.  The  number  density  slightly  increases  in 
a  layer  centered  around  20  km  because  of  the  thermal  and 
dynamic  nature  of  the  atmosphere  [10].  Tropospheric  aerosol 
particles  play  an  important  role  in  the  attenuation  of  solar 
radiation  because  they  strongly  affect  the  Earth’s  radiation  bud¬ 
get.  The  influence  of  aerosols  on  radiation  that  passes  through  the 
atmosphere  cannot  be  neglected,  especially  in  urban  or  indus¬ 
trialized  areas,  [20,21,23]. 

Aerosols  are  a  minor  constituent  of  the  Earth’s  atmosphere, 
but  they  play  an  important  role  in  the  energy  balance  of  the 
earth-atmosphere  system.  Aerosols  can  influence  the  climate  in 
two  ways  [24].  The  magnitude  of  these  influences,  however,  is 
poorly  constrained  because  of  limited  knowledge  on  the  pro¬ 
cesses  that  control  the  distributions  as  well  as  the  physical, 
chemical,  and  optical  properties  of  aerosols  [25],  First,  aerosol 
particles  can  reduce  solar  energy  by  absorbing  or  scattering  the 
radiation  under  cloud-free  and  cloudy  conditions,  thereby  affect¬ 
ing  the  energy  balance  of  the  Earth  [24].  This  direct  effect  is 
related  to  the  optical  properties  of  aerosol  particles,  such  as  single 
scattering  albedo,  size  distribution,  size  variability  with  relative 
humidity,  complex  refractive  index,  and  solubility.  The  levels  of 
scattering  and  absorption  depend  on  the  physical  and  chemical 
characteristics  of  aerosol  particles.  Consequently,  aerosols  modify 
the  Earth’s  radiation  budget  and  consequently  influence  the 
warming  or  cooling  of  the  planet.  Second,  aerosols  alter  the 
formation  and  the  precipitation  efficiency  of  liquid-water,  ice, 
and  mixed-phase  clouds,  thereby  causing  indirect  radiative  for¬ 
cing  associated  with  these  changes  in  the  cloud  properties.  The 
quantification  of  indirect  radiative  forcing  caused  by  aerosols  is 
an  especially  difficult  task  [26]. 

Radiative  forcing  of  the  climate  by  aerosols  may  be  compar¬ 
able  in  magnitude,  but  opposite  in  sign,  to  that  of  greenhouse 
gases.  Greenhouse  gases  trap  radiation  and  warm  the  surface  of 
the  earth,  while  aerosols  absorb  and  scatter  radiation  back  to 
space,  theoretically  resulting  in  a  net  cooling  effect  [24],  This 
variability  is  largely  due  to  the  much  shorter  atmospheric  lifetime 
of  aerosols  compared  with  that  of  greenhouse  gases  [23].  Turbid¬ 
ity  observations  had  been  made  since  1923,  but  their  measure¬ 
ments  are  inherently  difficult  and  expensive  even  making  some  of 
them  not  possible  because  of  the  strong  influence  of  other  atmo¬ 
spheric  components  [27,15,28].  The  accurate  determination  of 
turbidity  normally  requires  clear-sky  spectral  radiation  data 
obtained  using  sun  photometers  or  spectroradiometers.  These 
instruments  are  expensive  and  scarce,  and  thus,  turbidity  is 
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instead  generally  estimated  from  broadband  irradiance  measure¬ 
ments  [29]. 


3.  Atmospheric  turbidity  indices 

Aerosols  have  a  relationship  with  the  attenuation  of  solar 
radiation  that  reaches  the  Earth’s  surface,  and  thus,  a  number  of 
methods  and  atmospheric  turbidity  indices  have  been  introduced 
and  developed  in  the  past  decades  to  evaluate  atmospheric 
turbidity  and  to  quantify  its  influence  on  the  direct  irradiance 
on  the  Earth's  surface  [14,30,31].  The  Angstrom  turbidity  coeffi¬ 
cients  a  and  p  and  the  Schiiepp  turbidity  coefficient  B,  among 
others,  are  spectral  in  nature,  such  that  they  reference  the 
spectral  aerosol  optical  depth  to  a  common  wavelength  (i.e„ 
1000  nm  for  p  and  500  nm  for  B ).  On  the  other  hand,  the 
Unsworth-Monteith  turbidity  coefficient  Tu  and  the  Link  turbidity 
factor  Tl  are  broadband  in  nature,  i.e„  they  are  estimated  by 
reversing  or  inverting  the  models  that  estimate  broadband  beam 
irradiance  upon  measurement.  The  Unsworth-Monteith  turbidity 
coefficient  is,  by  definition,  equal  to  the  broadband  aerosol  optical 
depth,  which  is  the  sum  of  aerosol  absorption  and  scattering 
coefficients  integrated  over  the  entire  depth  of  the  atmosphere 
[15,30],  The  definition  of  each  turbidity  index  is  explained  in  the 
succeeding  sections. 

3.3.  Angstrom’s  turbidity  coefficients  (a  and  P) 

Under  clear  sky  conditions,  atmospheric  aerosols  generally 
contribute  to  the  largest  attenuation  of  solar  radiation  in  the 
visible  range.  Scattering  and  absorption  effects  cause  the  attenua¬ 
tion.  Since  attenuation  effects  of  the  scattering  and  absorption  by 
aerosols  are  difficult  to  separate,  Angstrom  suggested  a  single 
formula  that  is  generally  known  as  the  Angstrom  turbidity 
formula,  which  can  be  expressed  as  follows  [10,21,32]: 

=  (1) 

where  xa2  is  the  spectral  aerosol  optical  depth  at  the  wavelength, 
P  is  the  Angstrom  turbidity  coefficient,  a  is  the  wavelength 
exponent,  and  /  is  the  wavelength  in  micrometer. 

In  Eq.  (1),  p,  which  varies  from  0.0  to  0.5  or  even  higher,  is  an 
index  that  represents  the  amount  of  aerosols  in  the  atmosphere  at 
the  vertical  direction.  The  wavelength  exponent  a  is  related  to  the 
size  distribution  of  the  aerosol  particles.  Large  values  of  a  indicate 
a  relatively  high  ratio  of  small  to  large  particles,  which  vary  from 
0  to  4.  When  the  aerosol  particles  are  very  small  in  relation  to  air 
molecules,  a  should  approach  4,  whereas  for  very  large  particles, 
a  should  approach  0.  Generally,  a  is  has  a  value  between  0.5  and 
2.5.  The  average  a  for  the  entire  solar  spectrum  has  a  value 
between  0.9  and  2.0.  An  acceptable  average  value  for  most  natural 
atmospheres  is  a=  1.3 +  0.5.  Table  1  lists  the  combinations  of 
parameters  p  and  a,  which  may  be  considered  for  various  degrees 
of  atmospheric  cleanliness  [10,33], 

The  Angstrom  exponent  a  is  a  widely  used  parameter  in 
atmospheric  sciences  that  examines  the  optical  properties  of 
aerosol  particles.  Since  the  early  [34,35]  and  later  [36,37]  pub¬ 
lications  of  Angstrom,  where  the  a  parameter  was  mainly  used  to 


Table  1 

Parameters  for  various  degrees  of  atmospheric  cleanliness  [10], 


Atmosphere 

p 

a 

Visibility  (km) 

Clean 

0.00 

1.30 

340 

Clear 

0.10 

1.30 

28 

Turbid 

0.20 

1.30 

11 

Very  turbid 

0.40 

1.30 

<5 

describe  the  spectral  behavior  of  atmospheric  extinction  and 
transmission,  a  is  now  also  applied  to  a  variety  of  similar  but 
slightly  different  optical  properties  such  as  the  atmosphere 
scattering  or  backscattering  coefficients.  The  Angstrom  exponent 
is  very  popular  because  of  the  simplicity  of  its  equation  that 
enables  the  interpolation  or  extrapolation  of  aerosol  optical 
properties  and  because  of  its  connection  to  particle  microphysics 
(related  to  the  mean  size  of  aerosols)  [38]. 

3.3.3.  Techniques  for  determining  the  Angstrom  turbidity  coefficients 

(ct  and  p) 

Angstrom  turbidity  coefficients  a  and  p  can  be  determined 
using  a  number  of  techniques.  Traditionally,  a  and  p  are  derived 
from  either  spectral  direct  solar  radiation  measurements  or  from 
broadband  direct  solar  radiation  and  precipitable  water  measure¬ 
ments  [39,40]. 


3.3. 3. 3.  The  Volz  method.  Volz  [41]  applies  Eq.  (1)  to  sun  photometric 
measurements  performed  in  two  different  narrow  bands  that  are 
chosen  to  minimize  the  influence  of  other  atmospheric  attenuators, 
such  as  water  vapor.  In  this  case,  a  can  be  obtained  as  follows: 


g  dlm«A  lnG“0 

d  ln2  ~  in£ 


(2) 


where  X^  and  X2  represent  the  effective  wavelengths  in  the 
corresponding  narrow  bands  [21,26]. 

The  chosen  wavelengths  are  usually  380  and  500  nm,  where 
weather  conditions  play  a  very  important  role  in  the  distribution 
of  pollutants  in  the  air  [39,28].  At  380  nm,  no  molecular  absorption 
is  observed,  while  at  500  nm,  the  ozone  has  weak  absorption  [10]. 
Eq.  (2)  indicates  that  a  is  the  negative  of  the  slope  or  the  negative  of 
the  first  derivative  of  xx2  versus  the  wavelength  in  logarithmic  scale. 
After  a  is  derived,  p  can  be  calculated  in  the  same  wavelength 
interval  using  the  values  of  a  and  xa2  at  a  certain  X  [21].  However,  p 
alone  can  be  measured  with  a  single-wavelength  Volz  instrument  by 
assuming  that  a=1.3.  Alternately,  p  can  be  measured  at  2=1  pm, 
where  the  exponent  a  is  not  included  in  the  calculations.  At  this 
wavelength,  only  the  water  vapor  has  a  weak  absorption  band.  This 
method  is  very  accurate  and  yields  the  values  of  a  and  p  simulta¬ 
neously  [10,39,28]. 


3. 3. 3.2.  Direct  method.  The  direct  method  enables  the  calculation 
of  the  Angstrom  turbidity  coefficients  a  and  p  by  using  the 
following  formula  obtained  from  Eq.  (1)  [42]: 


d[ Inx^}  X  [ dx  ctk 

dln2  —  xa2  dX 


(3) 


The  parameter  p  can  also  be  determined  from  Eq.  (1). 

3. 3. 3.3.  Linear  fitting.  Linear  fitting  depends  on  the  omission  of  xai 
in  the  strong  absorption  bands  of  water  vapor  and  mixed  gases, 
and  it  determines  the  Angstrom  parameters  at  remaining 
wavelengths.  This  fact  constitutes  the  main  difference  with  the 
previous  (direct)  method.  The  Angstrom  parameters  are  derived 
by  applying  the  least-squares  fit,  as  written  in  Eq.  (1),  to  the 
remaining  wavelengths,  which  can  be  rewritten  in  the  following 
linear  form  [42]: 

IriXaA  =  -a  lnA+\np  (4) 

In  this  approximation,  the  slope  of  the  straight  line  yields  a, 
while  the  intercept  provides  /?.  Hence,  a  is  the  negative  of  the 
slope  (first  derivative)  of  In  xa2  versus  In  X.  This  method  has  often 
been  used  for  determining  aerosol  characteristics  because  of  its 
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simplicity  [21,43].  Since  the  Angstrom  formula  does  not  exactly 
fit  the  data  points,  a  second-order  polynomial  fit  for  In  Ta2  versus 
In  X  data  provides  excellent  agreement  with  the  measured  Ta2, 
resulting  in  order  differences  of  the  uncertainty  of  the  measure¬ 
ments  (~  0.01 -0.02),  while  a  linear  fit  yields  significant  differ¬ 
ences  with  the  measured  Ta2.  The  second-order  polynomial  fit  is 
applied  to  the  ta/l  values  at  all  wavelengths  and  is  expressed  in 
the  following  form  [26]: 

lnT^  =  a2  \nX2  +  oi]  In2+a0  (5) 

Since  a  is  the  negative  of  the  first  derivative  of  In  versus  In  X, 
the  following  expression  can  be  derived  from  Eq.  (4)  to  quantify 
the  spectral  dependence  of  or. 

a  = -2a2  ln/l-a!  (6) 


3. 1.1.4.  Horizontal  visibility.  Routine  measurements  of  p  and  a 
with  a  dual-wavelength  sun  photometer  are  conducted  at  a 
number  of  locations.  When  such  measurements  are  unavailable, 
the  turbidity  parameter  p  may  be  determined  from  the 
measurement  of  visibility  (Vis)  in  the  horizontal  direction,  if 
such  a  measurement  exists.  This  visibility  is  also  called  the 
meteorological  range.  For  visibilities  greater  than  5  km,  the 
turbidity  parameter  can  be  determined  using  the  following 
equation  developed  by  McClatchey,  Selby,  and  others  [10]: 

£  =  0.55*(3.912/Vis— 0.01 162)[0.02472(Vis—5)  +  l. 132]  (7) 

where  Vis  is  in  kilometers.  In  order  to  predict  /?,  the  value  of  a 
must  first  be  known.  Unfortunately,  even  with  the  help  of  proper 
instruments,  the  measured  visibility  is  slightly  subjective  in 
nature.  Therefore,  the  relationship  above  only  yields  a  gross 
estimate  of  p. 

3. 3. 1.5.  Rangarajan  and  Mani  method.  Rangarajan  and  Mani  [44] 
described  a  new,  simpler,  and  less  expensive  method  of  evaluating 
/ 3  from  direct  and  diffused  solar  radiation  measurements.  Their 
method  has  been  a  routine  in  principal  radiation  stations. 


3d  0.5Sa  +  0.75Sd 
lb  ”  (1— Sa)(l— Sd) 

(8) 

sd=\-g(pr 

(9) 

0.606m* 

(10) 

S““  643 +  m* 

g(/?)  =  1  —  1.375/3  for  p<  0.08 

(ID 

g(P)  =  0.9843-1.1429^  for  p>  0.08 

(12) 

where  m*  is  an  absolute  air  mass,  lb  is  the  direct  solar  radiation, 
and  Id  is  the  diffused  solar  radiation. 


C  =  1.003— 0.125a 
D=  1.089a +  0.5123 

where  E0  is  a  correction  term  for  the  earth-sun  distance. 


3. 3. 3.7.  Gueymard  model.  Gueymard  [47,48]  developed  a  new 
model  for  solar  beam  spectral  irradiance.  This  model  uses  the 
following  two  possible  a  values:  oq  for  X  <  0.5  pm  and  a2  for  other 
measurements.  The  p  parameter  is  independent  of  the 
wavelength  and  is  derived  from  the  following  relation: 

fie  =  0-5a2  |t5  (14) 

where  t5  is  the  aerosol  optical  thickness  at  0.5  pm.  Consequently, 
only  a  single  measurement  of  the  beam  spectral  irradiance  at  this 
wavelength  is  necessary. 


3. 3. 3.8.  Pinazo  et  al.  method.  Pinazo  et  al.  [40]  presented  a  new 
method  (referred  to  as  Method  II)  for  the  calculation  of  the 
Angstrom  turbidity  coefficient  based  on  the  ratio  of  direct  solar 
radiation  to  global  solar  radiation  on  a  horizontal  surface  and  on 
the  “C”  model  of  Iqbal.  Method  II  permits  the  calculation  of  ft 
without  employing  certain  variables,  such  as  the  thicknesses  of 
the  precipitable  water  and  of  the  ozone  layer.  Furthermore, 
measurements  on  dry  temperature  and  relative  humidity  are 
not  necessary.  According  to  Method  II,  ft  can  be  expressed  as 
follows: 


fipin  —  Itl 


(D— 0.1 2445a +  0.01 62)/(l. 003— 0.125a) 
mr(1.089a  +  0.5123) 


(15) 


The  definition  of  D  can  be  found  in  the  original  paper  of  the 
author  [40]. 


3. 3. 3.9.  Gueymard  and  Vignola  method.  Gueymard  and  Vignola 
[29]  proposed  a  semi-physical  method  to  evaluate  turbidity 
from  broadband  irradiance  measurements  and  from  other 
atmospheric  parameters.  This  method  should  be  under  cloudless 
skies  and  requires  very  accurate  measurements  of  global  (or 
diffused)  and  direct  radiations.  From  the  number  of  parametric 
runs  performed  with  SMARTS2,  I<db  {Kdb=ld!lb the  ratio  between 
diffused  and  direct  beam  normal  irradiance)  can  be  accurately 
parameterized  as  follows: 

Kdb  =  (a0+a,p+a2p2)/(L+a3p2')  (16) 

where  the  coefficients  a,  are  functions  of  the  zenith  angle, 
pressure  ozone  amount,  and  precipitable  water. 

If  an  experimental  value  of  I<db=ld/Ibn  is  obtained  from  the 
measurements,  Eq.  (16)  can  easily  be  solved  to  obtain  ft  as 
follows: 

P  =  0.5  j  [a’i-4(a2-a3Kdb)(a0-l<db))05~a-l  }/(a2-a33Cd6)  (17) 


3.3. 3.6.  Louche  et  al.  method.  Louche  et  al.  [39]  described  a  method  for 
determining  the  Angstrom's  turbidity  coefficient  in  a  Mediterranean 
location  from  the  total  pyrheliometric  measurements  in  the  para¬ 
meterization  of  the  direct  solar  total  irradiance  by  Bird  and  Hulstrom 
[45,46],  which  can  be  expressed  as  follows: 

B=  -^lnU^  (13) 

r  mrD  \A-BJ  y 

where  mr  is  the  relative  air  mass  at  actual  pressure  (dimensionless) 
and 

A  =  3ta/(0.975£03ScTrT0TgTw) 


B  =  0.1244a— 0.0162 


3.3.3.30.  Malik  method.  Malik  [7]  modified  the  Pinazo  et  al.  model 
[40]  to  introduce  a  new  variable,  called  the  “horizontal  visibility,” 
to  more  accurately  estimate  the  turbidity  coefficient  of  Angstrom 
P  from  the  ratio  of  horizontal  direct  solar  irradiance  to  the  global 
horizontal  irradiance  and  horizontal  visibility. 


P  =  -In 


(D— 0.1 2445a +  0.91 62)/(l. 003— 0.125a) 
mr(1.089a  +  0.5123) 


(18) 


3.3.3.)  3.  Utrillas  et  al.  method.  Utrillas  et  al.  [19]  developed  a  new 
method  aimed  at  determining  the  Angstrom  turbidity  coefficients 
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by  measuring  broadband  direct  irradiance  by  using  filters. 

ln%i  =  In/?- a  \nXaff  (19) 

where  t is  the  aerosol  optical  thickness  at  a  specified  spectral 
band  and  2a//  is  the  effective  wavelength  that  corresponds  to  each 
filter. 

3.2.  Linke’s  turbidity  factor  (TL) 

Linke  [58],  [59]  introduced  the  Linke  turbidity  factor  ( TL )  and 
showed  the  number  of  dry  atmospheres  necessary  to  cause  the 
attenuation  of  extra-terrestrial  radiation  produced  by  the  real 
atmosphere  [49].  TL  refers  to  the  entire  spectrum  or  to  the  overall 
spectrally  integrated  attenuation  [50].  Gueymard  [51]  suggested 
that  the  TL  coefficient  is  not  a  pure  turbidity  coefficient  because  it 
also  incorporates  water  vapor  and  N02  optical  depths.  Never¬ 
theless,  the  Tl  coefficient  may  be  useful  when  no  data  on 
precipitable  water  are  available,  making  the  determination  of  other 
coefficients  difficult  ([52]reported).  Although  TL  does  not  explicitly 
show  the  effects  of  aerosols  and  water  vapor,  its  easy  measurement 
made  it  quite  popular  in  meteorological  stations  [53]. 

Tl  is  a  key  input  for  several  models  that  assess  the  down- 
welling  irradiance  under  clear  skies.  These  models  are  used  by 
several  communities  in  the  fields  of  renewable  energies,  clima¬ 
tology,  agro-meteorology,  and  atmospheric  pollution  [54].  TL  is 
also  useful  for  the  prediction  of  the  availability  of  solar  radiation 
and  daylight  illuminance  under  cloudless  skies  [55].  Although  TL 
is  a  useful  parameter  in  many  applications,  it  has  a  serious 
drawback.  TL  varies  with  air  mass  even  when  atmospheric  condi¬ 
tions  remain  constant  [39].  TL  normally  varies  from  1  to  10  [31].  In 
a  pure  Rayleigh  atmosphere,  TL=1.  The  value  closest  to  this  ideal 
value  is  achieved  in  extremely  clear  and  cold  air  at  high  latitudes, 
Tl= 2.  In  a  polluted  atmosphere,  TL  can  increase  to  8  [56],  The 
attenuation  of  the  radiation  under  a  clear  sky  atmosphere 
becomes  larger  with  increasing  TL. 

3.2 A .  Techniques  for  determining  the  Linke's  turbidity  factor  (TL) 

Tl  can  be  directly  obtained  from  the  pyrheliometric  measure¬ 
ments  of  the  beam  irradiance  during  periods  of  veiy  clear  sky  with 
an  air  mass  of  2.  Since  this  kind  of  experimental  data  is  rarely 
available,  TL  generally  becomes  an  estimated  parameter.  In  addition, 
the  time  series  of  the  radiation  data  is  generally  too  short  to  allow 
estimation  on  a  daily  basis.  Fortunately,  long-term  monthly  average 
values  TLm  are  sufficient  for  most  applications.  Several  methods  for 
estimating  Tliti  values  can  be  found  in  literature  [27,57]. 

3.2.1A.  Linke.  Linke  [58]  defined  the  integral  Rayleigh  optical 
thickness  SR(mr )  as  the  integrated  optical  thickness  of  the 
terrestrial  atmosphere  that  is  free  of  clouds,  water  vapor,  and 
aerosols.  He  calculated  dR(mr)  from  theoretical  assumptions  and 
validated  it  in  a  very  pure  and  dry  mountain  atmosphere  by  using 
the  following  formulation: 

5R(mr)  =  0.128-0.054  logmr  (20) 


3.2. 3.2.  Feussner  and  Dubois.  Feussner  and  Dubois  [60]  published 
a  series  of  spectral  data  tables  that  enables  the  calculation  of 
SR(mr),  where  both  molecular  scattering  and  the  absorption  by 
the  stratospheric  ozone  layer  are  taken  into  account. 

3.2. 3. 3.  Valko.  Valko  [61]  suggested  empirical  fits  of  TL  as  a 
function  of  aerosol,  water  vapor  content,  and  relative  air  mass. 


Valko’s  formula  valid  for  5°  <  y  <  65°can  be  expressed  as  follows: 
Tl  =  (8  +  0.54) [1.75  log(w/mr  +  0.1)  +  14.5]-5.4  (21) 

where  8=1.07/3  (by  Schiiepp). 

3.2A.4.  Dogniaux.  Dogniaux  [62]  derived  the  following  relation 
from  extensive  experimental  campaigns,  illustrating  the  observed 
variation  in  TL  with  solar  height  y  (°),  atmospheric  water  vapor 
content  w  (cm),  and  the  Angstrom’s  turbidity  coefficient  /3: 

TL  =  (39  5g5w^47  4  +01)  +06  +  0.22 w)/3  for  5°  <  7  <  65“  (22) 

3.2A.5.  Jaenicke  and  Kasten.  Jaenicke  and  Kasten  [63]  developed 
the  following  linear  relationship  to  determine  the  atmospheric 
turbidity  factor  from  the  threshold  solar  elevation  that  is 
measured  by  any  given  sunshine  autograph  with  known  threshold 
irradiance  Ics: 

TL  =  0.154  ln(30/3cs)oc+1.05  ln(30/3cs)  for?  30°  (23) 

where  30  is  the  irradiance  of  the  extra-terrestrial  solar  radiation  and 
f0  =  1375  W  m  2  is  the  mean  sun-earth  distance. 

3.2A.6.  Kasten.  To  facilitate  the  determination  of  TL(mr)  from 
pyrheliometric  measurements,  Kasten  [64]  presented  the  following 
simple  analytic  formula  for  SR(mr),  known  as  Kasten’s  pyrhelio¬ 
metric  formula,  based  on  numerical  values  published  by  Feussner 
and  Dubois  [59]  and  on  international  recommendations  by  Special 
Committee  for  the  International  Geophysical  Year  (CSAGI  1958): 

SR(mr)  =  (9.4+0.9mr)_1  (24) 

3.2A.7.  Katz  et  al.  Katz  et  al.  [65]  experimentally  discovered  that 
Tl=2. 0+19/3,  but  they  used  the  definition  by  Kasten  [64]  on  the 
integrated  optical  thickness  of  the  terrestrial  atmosphere,  which 
only  includes  Rayleigh  scattering  and  ozone  absorption. 

3.2A.8.  Kasten.  TL  at  height  z  meters  above  sea  level  Tt(z)  was 
objectively  defined  by  Kasten  [66]  as  ([61]  reported)  follows: 

rt(z)=  1  +  (Sd(z)/<5r(z)  (25) 

where  <5 R(z)  is  the  relative  optical  thickness  associated  with 
Rayleigh  scattering  by  the  gaseous  molecules  in  the  atmosphere 
and  ozone  absorption,  and  dD(z)  is  the  relative  optical  thickness 
associated  with  aerosol  extinction  and  gaseous  absorption  other 
than  the  ozone  in  the  stratosphere. 

3.2A.9.  Louche  et  al.  Louche  et  al.  [50]  used  more  recent  data  on 
the  spectral  extra-terrestrial  solar  radiation  and  on  the  spectral 
absorption  coefficients  of  the  atmospheric  gases,  as  compiled  by 
Iqbal  [67],  to  calculate  the  improved  numerical  values  of  the 
integral  Rayleigh  optical  thickness  SR(mr)  according  to  the 
integration  procedure  outlined  by  Kasten  [64].  Based  on  the 
new  values  of  SR(mr),  Louche  et  al.  proposed  a  new  analytic 
formula  for  the  function  SR(mr )  and  demonstrated  that  the  new 
SR(mr )  values  are  definitely  higher  than  the  old  CSAGI  values  used 
by  Kasten,  particularly,  at  small  relative  optical  air  masses  mr. 
Louche  et  al.  also  proposed  an  approximate  equation  for 
converting  old  turbidity  factors  TL(mr)  into  new  ones. 

SR(mr)  =  (6.5567  +  1.7513mr-0.1202m^  +  0.0065m^-0.00013m^)_1 

(26) 

3.2AA0.  Abdelrahman  et  al.  Abdelrahman  et  al.  [53]  conducted 
turbidity  measurements  in  Dhahran,  Saudi  Arabia  which  has  a 
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semi-arid  climate.  The  inter-relationships  between  TL  and  f!  were 
found  to  be  linear,  similar  with  that  in  Avignon  (France)  and  in 
Potsdam  (West  Germany).  However,  the  regression  coefficients 
were  different  for  Dhahran. 

P  =  — 0.1326  +  0. 0433Tt  (27) 

3.2 A  A  I.  Grenier  et  al.  Grenier  et  al.  [68]  proposed  a  fourth-order 
polynomial  for  the  function  SR(mr)  that  is  similar  with  Louche 
et  al.  [50]  but  with  slightly  different  coefficients.  They  confined  the 
range  of  validity  to  1  <  mr  <  6  to  achieve  a  better  fit  in  this  range. 

SR(mr)  =  (5.4729 +  3.0312mr-0.6329m;:  +  0.0910m;?-0. 0051 2m?)  1 

(28) 


Furthermore,  Grenier  et  al.  illustrated  the  dependence  of  TL 
with  air  mass  and  proposed  a  Linke  turbidity  index  that  was 
normalized  at  an  air  mass  of  2,  which  they  called  TLAM2<  and  was 
therefore  independent  of  air  mass.  They  also  investigated  the 
dependence  of  Tlamj  on  water  vapor  and  on  the  Angstrom  spectral 
turbidity  factor  and  presented  these  dependencies  in  the  form  of 
various  tables  and  polynomial  fits. 


P  =  - 0. 120  + 0.07  1TMm2  0.5  <w 
P  =  -0.128  +  0.0707^2  1 .2  <  w 

P  =  — 0.1 38  +  0.0697^1^2  w  >  2.8 


<1.6  cm 

(29) 

<  3.2  cm 

(30) 

cm 

(31) 

3.2AA5.  Cucumo  et  al.  Cucumo  et  al.  [71]  developed  a  simple  and 
general  calculation  model  for  the  determination  of  the  average 
daily  Linke  turbidity  factor  in  places  where  the  daily  global 
irradiation  on  the  horizontal  plane  is  known.  Their  model  was 
validated  using  experimental  data  for  the  turbidity  factor  TL 
obtained  in  Arcavacata  di  Rende  (Cosenza)  and  in  Casaccia 
(Rome).  The  model  was  also  experimentally  validated  using 
experimental  Tt  data  obtained  in  Arcavacata  di  Rende  and  in 
Casaccia.  Moreover,  polynomial  correlations  of  the  turbidity 
factor  had  been  developed  for  28  Italian  localities  where  the 
European  Solar  Radiation  Atlas  (ESRA)  provided  maximum 
radiation  values  during  clear  days.  Using  these  correlations,  the 
average  turbidity  factor  on  each  day  of  the  year  was  calculated. 
The  suggested  formula  for  estimating  Tt  is  given  by 

T  -lnUtaHhK/scfio) 

Il=  - 5 -  (36) 

dRmr 

where  Ihn  is  the  direct  hourly  normal  energy  (i.e.,  the  incident 
hourly  energy  on  a  plane  perpendicular  to  the  Sun’s  rays),  dR  is 
the  result  of  the  formula  by  Louche  et  al.  [39]  corrected  by  Kasten 
[69],  and  mr  is  the  expression  by  Louche  et  al.  [39]. 

3.2AA6.  Geiger  et  al.  Geiger  et  al.  [72]  summarized  the  equations 
that  describe  the  modified  European  Solar  Radiation  Atlas  (ESRA) 
model,  which  was  initially  presented  by  Rigollier  et  al.  [73]  and  then 
modified  by  Remund  and  Page  [74]  to  take  into  account  the 
variation  in  the  integral  optical  depth  with  the  altitude  of  the  site. 


3.2AA2.  Molineaux  et  al.  Molineaux  et  al.  [39]  noted  that  the 
expressions  by  Louche  and  Grenier  for  SR(mr)  became  divergent 
for  air  mass  measurements  that  were  greater  than  20  and  7.  They 
adopted  the  coefficients  of  the  original  expression  by  Linke  to 
take  into  account  the  absorption  of  permanent  gases  as  follows: 

5R(mr)  =  0.124-0.0656  logmr  (32) 

3.2 A  A3.  Kasten.  To  compare  and  judge  the  two  parameterization 
formulae  for  SR(mr)  proposed  by  Louche  et  al.  [50]  and  by  Grenier 
et  al.  [68],  respectively,  Kasten  [69]  presented  a  new  computation  of 
SR(mr )  according  to  his  outlined  integration  procedure  described  in 
[64],  which  was  repeated  by  Louche  et  al.  This  new  computation 
disregarded  several  simplifications  and  assumptions  made  by 
Louche  et  al.  and  by  Grenier  et  al.,  while  additional  parameters 
were  considered  and  some  recently  improved  input  data  were  used. 

TLa(mr)  =  TLK(mr)[\/5Ra(mr)}/[l/5RK(mr)}  (33) 

where  TLa(mr )  is  the  new  Linke  turbidity  factor,  TLK(mr)  is  the  Linke 
turbidity  factor  determined  by  the  algorithm  of  Kasten  [64],  <5jydmr) 
is  the  integral  Rayleigh  optical  thickness  calculated  using  the 
formula  of  Kasten  [64],  and  SRa(mr)  is  the  integral  Rayleigh  optical 
thickness  calculated  using  the  formula  of  the  new  spectral 
integration,  which  is  given  by: 

SRa(mr)  =  (6. 6296+1. 7513mr-0.1202m;:  + 0.0065m]? -0. 00013m]4) 

(34) 


3.2AA7.  Ineichen  and  Perez.  Ineichen  and  Perez  [75]  presented  a 
new  formulation  of  turbidity  that  displays  a  relatively  stable 
behavior  during  a  specific  day.  The  new  turbidity  coefficient 
agreed  with  that  by  Kasten  [64]  at  an  air  mass  of  2.  However, 
the  contribution  of  absorption  by  atmospheric  gases  into  the 
integral  optical  depth  was  not  taken  into  account.  The 
formulation  of  TLIP  can  be  expressed  as  follows: 

Tup  =  1  +  [ll.lln(b£0/f,n)]/mr  (37) 

where  b  is  given  by 

b  =  0.664 +0.1 63/exp (station  height,  meters/8000)  (38) 

The  new  formulation  has  the  following  advantages:  (1)  solar 
altitude-independent  and  (2)  matches  the  original  Linke  turbidity 
factor  at  the  air  mass  of  2.  Therefore,  the  new  formulation 
remains  coherent  with  previous  studies  and  is  based  on  high 
quality  data  that  are  representative  of  widely  differing  geographic 
locations,  altitudes,  and  climates. 

3.2 A  A  8.  Polo  et  al.  Polo  et  al.  [32]  presented  a  methodology  for 
estimating  daily  Linke  turbidity  factor  under  clear  sky  conditions 
at  an  air  mass  of  2  from  the  global  horizontal  irradiance  data  at 
solar  noon  and  from  monthly  mean  values  of  the  Linke  turbidity 
factor.  The  method  utilized  the  ESRA  clear  sky  model  to  establish 
a  functional  relationship  between  the  global  irradiance  at  solar 
noon  and  the  final  estimation  of  the  Linke  turbidity  factor. 


3.2.1  A4.  Kambezidis  et  al.  Kambezidis  et  al.  [70]  derived  an 
expression  for  the  Linke  turbidity  parameter  in  the  visible 
region  of  the  solar  spectrum  Tlvis ■  Their  expression  is  based  on 
newly  identified  visible  transmittances  for  various  atmospheric 
constituents  obtained  by  spectrally  integrating  wavelength- 
dependent  atmospheric  transmittances. 


ln((  bnvis )  +  Irk/sci/is'G) ) 
*5 Rvisttlr 


(35) 


3.3.  Schuepp's  turbidity  coefficient  (B) 


Schiiepp  [76]  introduced  the  turbidity  parameter  B,  which  is 
fundamentally  based  similarly  with  p,  but  referred  to  the  base  10 
instead  of  e  and  to  0.5  pm  wavelength  instead  of  1  pm. 


d,\/_  = 


log 


1  lo(Ai) 
Ep  /(At) 

mr 


— CEr  +  Tw) 


(39) 


where  /(Ax)  is  the  direct  solar  flux  on  the  Earth’s  surface,  /0(A/l)  is 
the  spectral  irradiance  above  the  atmosphere,  rR  is  the  mean 
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Rayleigh  extinction,  and  rw  is  the  selective  absorption  coefficient 
for  atmospheric  gases  02,  03,  H20,  C02,  and  so  on. 


3.3.  J.  Techniques  for  determining  Schuepp’s  turbidity  coefficient  (B) 
Katz  et  al.  [65]  identified  relationships  and  correlations  among 
various  atmospheric  turbidity  coefficients  based  on  the  direct 
radiation  measurements  in  five  bands  of  the  visible  region,  which 
were  performed  in  a  semi-rural  site.  Tables  2  and  3  present  the 
values  of  a  and  b  coefficients  in  the  linear  regression  Y=a+bX, 
which  were  calculated  among  different  turbidity  parameters  /?,  B,-, 
and  Bq.5- 


3.4.  Unsworth-Monteith's  turbidity  coefficient  (Tu) 


Unsworth  and  Monteith  [77]  introduced  a  new  method  for  the 
estimation  of  the  turbidity  of  the  atmosphere.  The  Unsworth- 
Monteith  turbidity  coefficient  Tu  represents  the  number  of  clean 
and  dry  atmospheres  (CDAs),  which  exhibit  the  same  attenuation 
as  the  actual  atmosphere  that  contains  aerosols  and  water  vapor. 
Typically,  Tu  varies  from  0  to  1.  Unsworth  and  Monteith  used 
experimental  irradiance  measurements  of  a  clean  atmosphere 
and  created  the  following  formulation: 

IbnA=ItnAexp(-Tmmr)  (40) 


where  Ibn  and  mr  are  defined  earlier,  Ibn  represents  the  normal 
incidence  total  solar  radiation  under  a  dust-free  atmosphere 
(W  m~2),  and  exp(-Tumr)  represents  an  aerosol  spectral  trans¬ 
mittance  coefficient.  With  the  definition  of  the  integrated  total 
optical  thickness,  an  integrated  expression  of  Tv  that  is  similar  to 
the  expression  of  TL  can  be  obtained  as  follows: 


Tu  =  ~ 


(41) 


3.5.  Illuminance's  turbidity  factor  (Tu) 


Illuminance  turbidity  factor  Tu  was  introduced  by  Navvab  et  al. 
The  concept  of  the  illuminance  turbidity  factor  is  analogous  to 
that  of  the  Linke  factor  with  Ibn,  I0,  and  TL  substituted  with  I vsn,  /„o, 


Table  2 

Values  of  a  in  the  linear  regression  Y=a+bX. 


Y/X 

P 

B, 

B2 

b3 

b4 

B5 

Bo.5 

P 

0 

-0.005 

0.008 

0.002 

-0.007 

-0.003 

-0.001 

B, 

0 

0.019 

0.018 

0.008 

0.010 

0.001 

b2 

0 

-0.003 

-0.010 

-0.009 

-0.004 

b3 

0 

-0.006 

-0.001 

0.006 

b4 

0 

0.002 

0.005 

b5 

0 

0.006 

Bq.5 

0 

Table  3 

Values  of  b  in  the  linear  regression  Y= 

a+bX. 

Y/X 

P 

B, 

b2 

b3 

b4 

b5 

Bo.5 

P 

1 

0.95 

1.39 

0.73 

0.87 

0.93 

1.02 

B, 

1 

1.36 

0.69 

0.83 

0.91 

1.07 

b2 

1 

0.51 

0.61 

0.68 

0.72 

b3 

1 

1.13 

1.24 

1.29 

b4 

1 

1.10 

1.19 

b5 

1 

1.05 

Bq.5 

1 

and  Tu-  Tu  is  apparently  a  better  parameter  for  describing  atmo¬ 
spheric  conditions  because,  unlike  the  Linke  factor,  it  is  insensi¬ 
tive  to  water  vapor  content,  air  mass,  and  solar  altitude  angle  or 
time  of  day.  Tu  can  be  determined  from  the  direct  measurement  of 
the  beam  normal  illuminance  by  Chaiwiwatworakul  et  al.  [31] 
and  Navvab  et  al.  [78].  Illuminance  turbidity  can  be  obtained  as 
follows: 


The  values  of  Su  can  be  calculated  as  follows: 
04 

"  (l  +  0.0045mr) 


(42) 


(43) 


Illuminance  turbidity  has  a  heuristic  advantage  as  an  approx¬ 
imate  analytical  expression  because  it  can  be  directly  derived 
from  the  Angstrom  formula  as  follows: 

Tu  =  1+21.6/?  (44) 


3.6.  Volz’s  turbidity  factor  (Tv) 


Computations  of  Volz  turbidity  are  performed  according  to  the 
usual  Bouguer-Lambert-Beer  Law  that  expresses  the  measured 
intensity  at  wavelength  X  as  follows  [79]: 


Tv  =  — 


lnl6n  +  ln£o~ln/o 


(45) 


where  is  the  Rayleigh  scattering  coefficient  for  air  molecules  at 
wavelength  X,  and  t°?z  is  the  absorption  coefficient  for  ozone  at  X. 


3.7.  Kastrov’s  turbidity  coefficient  (TK) 


Adeyefa  et  al.  [80]  reported  that  Kastrov  [81  ]  obtained  the 
following  simple  formula  for  the  integral  solar  radiation  on  the 
surface  level: 


4  (I+TKm*) 

from  which  the  following  can  be  obtained: 


(46) 


(47) 


where  Tc  is  a  quantitative  characteristic  of  atmospheric  transparency. 


3.8.  Rational's  turbidity  factor  (Tr) 


Majumda  et  al.  [82]  defined  a  new  measure  of  total  atmo¬ 
spheric  turbidity,  termed  as  the  rational  turbidity  factor,  Tr,  to 
overcome  the  limitations  of  the  Linke  turbidity  factor  TL.  Tr  is 
composed  of  three  components,  namely,  (i)  pare  and  dry  air  (the 
basic  effect),  (ii)  precipitable  water  vapor  w,  and  (iii)  aerosol 
particles  in  the  atmosphere  (dust,  smoke,  and  haze).  The  formula¬ 
tion  of  Tr  is  expressed  as  follows: 


1  70.32491  — log/to\ 1/057 
r  mr\  0.072375  J 


(48) 


Majumda  et  al.  [82]  correlated  the  rational  turbidity  factor  Tr 
with  the  Angstrom-Schiiepp  turbidity  coefficient  B  for  different 
values  of  precipitable  water  w  at  unit  air  mass  to  eliminate  the 
effect  of  any  possible  virtual  variation  with  air  mass.  They 
demonstrated  that  the  effect  of  B  and  w  were  inseparable  because 
the  scattering  of  radiation  by  aerosol  particles  and  the  absorption 
of  the  same  by  water  vapor  coexist  within  the  absorption  bands  of 
water  vapor.  Thus,  the  two  effects  cannot  be  isolated  from  each 
other.  This  interaction  becomes  more  pronounced  with  increasing 
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particle  size  since  the  negative  exponent  of  A  in  the  Angstrom- 
Schiiepp  equation  approaches  zero  for  large  particles. 

Majumda  et  al.  [83]  also  investigated  the  effects  of  varying  air 
masses  mr  on  Tr  in  relation  to  B  and  w  and  presented  a  formula 
that  connects  T,  with  B,  w,  and  mr  as  follows: 


observed  when  comparing  data  from  different  parts  of  a  day  or 
year. 


References 


Tr=  l+80(e/1.8)1+0S17<0586)mr  + 


25 


1-0.2347(0.869 )mr~ 


(49) 


4.  Conclusion 

The  results  of  this  review  indicate  the  following: 

1.  In  the  last  century,  a  number  of  atmospheric  turbidity  indices 
were  introduced  and  several  techniques  were  developed  to 
determine  their  values.  The  most  currently  used  indices 
include  the  Angstrom  turbidity  coefficient  p,  the  Linke  turbid¬ 
ity  factor  Tl,  the  Schuepp  turbidity  coefficient  B,  and  the 
Unsworth-Monteith  Tu. 

2.  The  Angstrom  exponent  a  is  very  popular  because  the  simpli¬ 
city  of  its  equation  enables  the  interpolation  or  extrapolation 
of  aerosol  optical  properties  and  because  of  its  connection  to 
particle  microphysics.  In  addition,  the  Angstrom  exponent  a 
also  indicates  the  amount  of  aerosols  and  represents  the 
combined  effects  of  scattering  and  absorption  caused  by 
aerosols. 

3.  The  Angstrom  turbidity  coefficient  ji  has  long  been  used  as  a 
turbidity  index  and  it  can  be  determined  by  different  methods 
for  spectral  and  broadband  radiation  measurements  in  large 
spectrum  bands  that  are  not  affected  by  water  vapor  absorp¬ 
tion  (280  nm  <  A  >  630  nm). 

4.  p  cannot  be  evaluated  unless  a  is  known  a  priori.  This 
uncertainty  in  a  makes  the  determination  of  p  less  precise 
and  discriminating. 

5.  The  Linke  turbidity  factor  TL  is  important  for  the  calculation  of 
the  beam  solar  and  the  diffused  radiation  in  a  specific  locality. 
This  parameter  is  also  useful  for  the  prediction  of  available 
solar  radiation  and  daylight  illuminance  under  cloudless  skies. 
The  Linke  turbidity  factor  TL  is  also  a  key  input  to  several 
models  that  assess  down-welling  irradiance  under  clear  skies 
and  that  are  used  by  several  communities  in  the  fields  of 
renewable  energies,  climatology,  agro-meteorology,  and  atmo¬ 
spheric  pollution. 

6.  The  main  problem  with  TL  is  that  it  is  difficult  to  be  accurately 
determined  and  it  fails  to  distinguish  the  effect  of  water  vapor 
from  the  absorption  caused  by  aerosol.  TL  is  also  not  truly 
independent  of  the  optical  air  mass  and  it  varies  diurnally  even 
under  an  unchanging  atmosphere.  Moreover,  TL  is  subject  to 
virtual  daily  variations,  such  as  broadband  aerosol  optical 
depth  5a  caused  by  the  combination  of  different  parasitic 
effects.  Virtual  variations  in  TL  caused  by  the  solar  zenith  angle 
y  are  smoother  compared  with  those  of  Sa  for  y  <  80°,  but 
become  steeper  for  larger  zenith  angles.  A  reference  air  mass 
value  of  2  should  be  calculated  to  remove  a  part  of  the  TL’s 
parasitic  variation. 

7.  The  Angstrom  and  the  Schuepp  turbidity  coefficients,  p  and  B 
respectively,  have  the  same  spectral  definition.  They  both 
correspond  to  at  1000  nm  (base  e)  and  at  500  nm  (base 
10),  respectively,  and  thus,  they  are  true  turbidity  coefficients 
that  are  unaffected  by  anything  except  by  the  aerosol  total 
burden. 

8.  The  Unsworth-Monteith  coefficient  (or  broadband  aerosol 
optical  depth)  was  found  to  be  slightly  dependent  on  both 
zenith  angle  and  water  vapor.  Hence,  caution  should  be 
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